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Abstract 

We introduce a generalization of Glimm's random choice method, which provides us with 
an approximation of entropy solutions to quasilinear hyperbolic system of balance laws. 
The flux-function and the source term of the equations may depend on the unknown as well 
as on the time and space variables. The method is based on local approximate solutions 
of the generalized Riemann problem, which form building blocks in our scheme and allow 
us to take into account naturally the effects of the flux and source terms. To establish the 
nonlinear stability of these approximations, we investigate nonlinear interactions between 
generalized wave patterns. This analysis leads us to a global existence result for quasilinear 
hyperbolic systems with source-term, and applies, for instance, to the compressible Euler 
equations in general geometries and to hyperbolic systems posed on a Lorentzian manifold. 

1. Introduction 

1.1 Hyperbolic systems of balance laws. This papeiH is concerned with the approxi- 
mation of entropy solutions to the Cauchy problem for a quasilinear hyperbolic system 

dtu + dxf{t,x,u) = g{t,x,u), t > 0, x G M, (1.1) 

u{0,x) =uoix), xgR, (1.2) 

where u = u{t, x) £ MP is the unknown. We propose here a generalized version of the Glimm 
scheme [TU] which allows us to deal with a large class of mappings /, g and take into account the 
geometric effect of the flux and source terms. Our scheme is based on an approximate solver for 
the generalized Riemann problem, based on an asymptotic expansion introduced by LeFloch 
and Raviart [17j . The approach provides high accuracy and stability, under mild restrictions 
on the equation and the data. 

In (|l.ip . the flux / = f{t,x,u) € W and the source-term g ~ g{t,x,u) S W are given 
smooth maps defined for all (t, x, u) € M+ x M. x U, where U is a small neighborhood of the 
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origin in MP, and the initial data uq : M U is a function with bounded total variation. We 
assume that the Jacobian matrix A{t, x, u) := ^{t, x, u) admits p real and distinct eigenvalues, 

Xi{t, X, u) < X2{t, x,u) < . . . < \p{t, X, u), 

and therefore a basis of right-eigenvectors rj{t,x,u) (1 < j < p), Finally, we assume that 
each characteristic field is either genuinely nonlinear (VAj(t, a;, u) • ri{t,x,u) ^ 0) or linearly 
degenerate {\7Xj{t,x,u) ■ ri{t,x,u) — 0). 

One important motivation for considering general balance laws (ll.ip comes from the theory 
of general relativity. In this context, the vector u typically consists of fluid variables as well as 
(first order derivatives) of the coefficients of an unknown, Lorentzian metric tensor. (See [S] [5] 
and the reference therein.) One can also freeze the metric coefficients and concentrate on the 
dynamics of the fluid. For instance, the compressible Euler equations describing the dynamics 
of a gas flow in general geometry read: 

OtP + Oxipv) = pv p, 

a a 

d,{pv) + d.Apv^ +p)^ -^ipv^) - ^ pv, (1.3) 

a a 

dt{pE) + d.,{pvE + pv) ^ ^—{pvE + pv)~^ pE 

a a 

where a = a{t,x) > can be regarded as the cross section of a time-dependent (moving) duct, 
and p, V, p{p,e), e, and E = e + are the density, velocity, pressure, internal energy, and 
total energy of the gas, respectively. The system (|1.3|) describes a situation where the fluid does 
not affect the variation of the duct; i.e. the function a{t,x) is given and, for simplicity, smooth. 
The system (|1.3I) is of the form (|I.ip with u = {p, pv, pE^ , f — f{u) — [pv, pv^ +p, pvE+pv)'^ 
and g — g{t, x, u) — —^^gi{u) — where gi{u) — {pv, pv'^, pvE + pv)'^ . 

We are interested in solutions to (|l.ip - (ll.2p which have bounded total variation in space 
for all times and satisfy the equations in the sense of distributions, together with an entropy 
condition [13 IB [TB] . In the special case that 

f^fiu), 5 = 0, 

the existence of global entropy solutions was established by Glimm [10], assuming that the 
initial data uo{x) has sufficiently small total variation. Recall that two main ingredients in 
Glimm's random choice method are (1) the solutions of Ricmann problems and (2) a projection 
step based on a sequence of randomly chosen points. 

Let us first indicate some of the earlier work on the subject. The system (II. ip with 

f = f{x,u), g^g{x,u). 



was treated in pioneering work by Liu |20[ 121] . via a suitable extension of the Glimm method: 
the approximate solutions are defined by pasting together steady state solutions, i.e., solutions 
V = v{x) of the ordinary differential equation 

^(/(x,i;)) = g{x,v). 

He established the existence of solutions defined in a finite interval of time [0, T) as long as 
either T or the norms of g and dg/du are sufficiently small. Next, assuming in addition 
that the eigenvalues of the matrix A{x,u) never vanish (so that no resonance takes place), Liu 
deduced a global existence result (with T = +oo). Steady-state solutions were also used in the 
work by Glimm, Marshall, and Plohr JJJ. 
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For more general mappings /, g, the existence for (ll.ip - (|1.2p is established by Dafermos and 
Hsiao [7] and Dafermos [SI [9]. They assume that fx{u*,t,x) = g{u*,t,x) = at some (equi- 
librium) constant state u* , hence u* is a solution of (|l.ip around which (|l.ip can be formally 
linearized. They also require that the linearized system satisfies a dissipative property. Their 
main result concerns the consistency and stability of a generalization of the Glimm method, 
yielding therefore the global existence of entropy solutions to (|l.ip . In [7j, the approximate 
solutions to the Cauchy problem on each time step are based on classical Riemann solutions 
with initial data suitably modified by both the source term g and the map 9 :— fx- 

Next, Amadori et al. [TJ [2] developed further techniques to establish the existence of solu- 
tions for a large class of systems having / = f{u) and g = g{x, u), and discussed Dafermos- Hsiao 
dissipative condition. For some particular systems (of two or three equations) the condition 
that the total variation be small can be relaxed; see for instance Luskin and Temple [22,, Groah 
and Temple pj], Barnes, LeFloch, Schmidt, and Stewart [3], and the references cited therein. In 
these papers, the decreasing of a total variation functional (measured with respect to a suitable 
chosen coordinate) was the key to establish the stability of the scheme. 

1.2 A new version of the Glimm method. In the present paper we provide an alterna- 
tive approach to Dafermos-Hsiao's method, and introduce a generalized version of the Glimm 
scheme for general mappings f,g. Integrability assumptions will be required (and discussed 
later on) on the matrix A and the mapping q : M_(. xM. xU ^ W defined by 

df 

q{t, X, u) := g{t, x, u) - —{t, x, u). (1.4) 

It should be emphasized that only this combination of the source and the flux will be important 
in our approach, which can be summarized as follows. 

First, we study the generalized Riemann problem associated with the system (jl.ip . i.e. the 
Cauchy problem with piecewise constant initial data. The existence of solutions defined locally 
in spacetime in a neighborhood of the initial discontinuity was studied in Li and Yu T8] and 
Harabetian [13| . Contrary to the case where f,g only depend upon the unknown u, no closed 
formula is available for the solutions of the generalized Riemann problem. We propose here an 
approximate Riemann solver, inspired by a technique of asymptotic expansion introduced by 
Ben-Artzi and Falcovitz [4] (for the gas dynamics equations) and LeFloch and Raviart jl7j (for 
general hyperbolic systems of balance laws) ; see also iO^ . 

Our scheme for solving approximately the generalized Riemann problem can be re-interpreted 
as a splitting algorithm (the hyperbolic operator and the source term being decoupled). Since 
an approximate (rather than an exact) solution to the generalized Riemann problem is used, 
it is crucial to establish an error estimate which we achieve in Proposition 2.1 below, under a 
mild assumption on the data uq, /, g. This estimate will be necessary to ensure the consistency 
of our generalized Glimm method. 

Second, we study the nonlinear interaction of waves issuing from two generalized Riemann 
problems, and establish a suitable extension of Glimm's estimates [lOj to the general system 
(jl.ip : cf. Proposition 3.3. This is a key, technical part of our analysis. 

Third, we introduce our scheme and prove its stability in total variation, under the assump- 
tion that the initial data uq has sufficiently small total variation and that the total amplification 
due to (the derivatives of) /, g to the total variation of the solution is sufficiently small; cf. The- 
orem 4.3. More precisely, we impose that 

d^A d^A dq 
dtdu dxdu du 

are sufficiently small in L-'^(M+ x M). 

Finally, we conclude with the convergence of the proposed scheme (Cf. Theorem 5.1) which 
yields the global existence of entropy solutions for the Cauchy problem (|l.ip - (|1.2p . The solution 
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satisfies an entropy inequality and has bounded total variation in x for all t > 0. Our results 
cover in particular the case 

dtu + d,f{u)^git), (1.5) 
for which global existence of entropy solutions is established under the sole assumption 

r + co 

/ \g{t)\dt«l. (1.6) 
Jo 

Without further restriction on the flux /, this condition is clearly necessary in order to apply 
the Glimm method, since, for instance in the trivial case p = 1 and / = 0, (jl.Sp reduces to the 
differential equation 

dtu = g{t). (1.7) 

On one hand, the condition p.6p holds if and only if every solution of (|1.7p remains close to 
a constant state, which is a necessary condition in order to apply the Glimm method. On the 
other hand, when one of the eigenvalues of the system (|l.ip vanishes, the amplitude of solutions 
could become arbitrarily large and the solutions would not remain bounded — except when the 
source term satisfies a "damping" property in time. 

As a direct application, the global existence of entropy solutions to (|1.3p follows, if the source 
g and its derivative ^ are sufficiently small in L^(R+ x R), which is the case, for instance, if 
the support of [at^ax) is sufficiently small. 

2. An approximate solver for the generalized Riemann problem 



In the present section, we introduce an approximate solution to the generalized Riemann 
problem associated with the system (|l.ll) . and we derive an error estimates (see Proposition 2.1 
below) . 

Given > 0, xq G K, and two constant states ul.ub, e M^, we consider the generalized 
Riemann problem, denoted by Rq{u]^,uii] Iq, xq), and consisting of the following equations and 
initial conditions: 

dtU + dxfit,x,u) — g{t,x,u), t>to,xeR, (2.1) 



u{0,x) = ( ^<^0' (2.2) 

^ ' [ Ur, X > Xo- 

Replacing / and g in (|2.ip by /(^o, xq, u) and 0, respectively, the problem Rg(ul,ur; to, Xq) 
reduces to the classical Riemann proWem, which we denote by Rc{uL,UR;to,Xo), that is the 
equations 

dtu + dxf{tQ,xo,u) ^0, u{t,x)eW, t>tQ, xeR (2.3) 

together with the initial data (2.2). This problem was solved by Lax under the assumption that 
the initial jump {ur — ul| be sufficiently small: the solution to Rc{ul, ur; Iq, xq) is self-similar 
(i.e. depends only on jEjf) ^-i^d consists of at most {p+ 1) constant states — uq, ui, . . . ,Up — 
Ur, separated by rarefaction waves, shock waves or contact discontinuities; see Figure 2.1. 

The following terminology and notation will be used throughout this paper. Let Wc — 
Wci£,', ul, Ur; to, Xq) be the solution of Rc{ul,ur; to, xq) with ^ = {x — xo)/{t — to). We say 
that the problem Rc{uL,UR;to,xo) is solved by the elementary waves (ui_i,Ui) {i — 1,. .. ,p) 
if each constant state Ui belongs to the i-wave curve Wi(ui_i) issued from the state Ui^i in the 
phase space, and (ui_i, Ui) is called an i — wave of Rc{ul, ur; to, xo). When the ^-characteristic 
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field is genuinely nonlinear, the curve yVi(ui_i) consists of two parts, the i-rarefaction curve 
and the j-shock curve issuing from if j-characteristic field is linearly degenerate, the curve 
yVi{ui-i) is a curve of i-contact discontinuities. Call the strength of the i-wave {ui^i,Ui) 
along the i-curve, so that, for a genuinely nonlinear z-field, we can assume that £i > if 
Ui) is a rarefaction wave, and < if Ui) is a shock wave. On the other hand, Si 

has no specific sign if (ui^i,Ui) is a contact discontinuity. 

Let £i{uL, ur; to,xo) denote the wave strength of the j-wave {ui-i, ui) in the Riemann prob- 
lem Rc{ul,ur] to,xo), and vector s = (ei, . . . , Ep) denote the wave strength of Rc{ul, uii;tQ,XQ) 
(so |e| is equivalent to the total variation of 14^c'(^; m/,, u^; to, 2^0))- In addition, we let a~ — 
Xi{ui-i,to, xq) and — Xi{ui,to, xq) be the lower and upper speeds of the i-rarefaction wave 
{ui-i,Ui) respectively, and ct^ be the speed of the z-shock or i-contact discontinuity. If the 
i-wave is a shock or a contact discontinuity we set = cr^ = cr^. 

From the implicit function theorem we deduce that the states Ui and the speeds af^ are 
smooth functions of ul, ur, to, and xq. Moreover, one can check that Ui — ul + 0{1)\ur ~ ul\ 
(z = 0, 1, . . . ,p), and, for an i-shock 

(Ti = Xt{ui^i;to,xo) + 0(1) \ut - Ui-i\, i = 1,2, . . . ,p, 

where 0(1) is bounded function possibly depending on ul, ur (z U, to > 0, and xq G R. 

Consider next the generalized Riemann problem on which a large literature is available 
[IHl [131 m El [IZ] ■ First, we recall [T^ that the solution of Rq{ul,ur; io, a;o) is piecewise smooth 
and has a local structure which is similar to the one of the associated classical Riemann problem 
Rc{ul, Ur; to,xo)- Following ^17j we consider an approximate Riemann solution of the problem 
Rg{ul,ur; to, Xq), denoted by Wait, x; ul,ur; to, xo) and defined by 

Wait, x; ul,ur; to, xo) = Wc{0 + {t - h) q{to, xo, Wc{0) (2-4) 
for t > to and a; G M. Here, the function q{t,x,u) is given by (|1.4p . and 

Wc{0 = Wc{^;uL,UR;to,xo). 

t — to 

Observe that the function WG{t,x;uL,UR]to,xo) is constructed as a superposition of the cor- 
responding classical Riemann solution Wc{£,',ul, UR;to,xo) and an asymptotic expansion term 
{t-to)q{to,xo,Wc{C}) (see Figure 2.2). 

Within a region where function Wc{C} is a constant, the function Wait, x; ul, ur; to, xq) is 
a linear function of t, namely, 

WG{t,x;uL,UR;to,xo)^Ui + {t-to)q{to,xo,Ui), a^^ < ^ < (t:[_^^ (2.5) 

for i = 0,1, ... ,p. By convention, aQ :— —00 and cTj^+i := +oo. Whenever there will be no 
ambiguity, we will use the notation Wait, x) or Wait, x; ul, ur) for Wait, x; ul, ur; to, xq). 

To describe the structure of WG{t, x; ul,ur; to, xo), it is convenient to say that the approx- 
imate solution WG{t,x,UL,UR;to,xo) consists of an i-wave {ui^i,Ui) if {ui^i,Ui) is an i-wave 
of the corresponding classical Riemann solution Wc{^; u^, ur; to, xq). 
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Figure 2.2 : Generalized Riemann solution (p=2), 
UL{t), ui{t), uji{t) are functions of t and u{t,£,) is constructed by (12. 4p . 

We now prove that the function Wait, x) defined in (|2.4p approximately solves the problem 
RciuLTUR'jtQ, xq), by evaluating the discrepancy between Wcitjx) and the exact solution of 
RG{uL,uii;to,xo). Given any s > and r > 0, and any function 9 : x R ^ R with 
compact support, we now show that the term 

A(s, r; 9) -.^ f f {Wo dt9 + f{t, x, Wq) dj + g{t, x, Wq) 9) dxdt (2.6) 

J to Jxo-r 

is of third order in r, s, provided that the condition (|2.7p holds. 



Proposition 2.1. Let 9 : M+ x M ^ M &e a compactly supported, function. Then, for 
every {to,xo) G R_|_ x R, UL,uii G Z//, and any positive numbers s,r satisfying the (Courant- 
Friedrichs-Levy -type) stability condition 

-sup|A,(<,a::,'a)| < 1 (2.7) 
r 

(the supremum being taken over 1 < i <p, {t,x) £ R4. x M, and u ^U), the 
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function Wait, x) — Wait, x;uL,UR]to, xq) satisfies 




A(s,r;0) = / WGito + s,-)eito + sr)dx- Waito, -^0, ■)(!■. 



\x 




+ / /(•,a;o + r,H/G(-,a;o + r))6'(-,a;o+r)dt 



(2.8) 



Jto 



/(•, xo - r, Wg{-, xq - r))0{-, xq - r)dt 



to 



+ 0{l){s^ + r^)is + r + \u„, - wl|)||0||ci, 



where A{s,r;d) is given in (12. 6p and \ \0\\ci — \ \0\\co + \ \dtO\\co + \ \dxd\\co . 

The left-hand side of (12. 8p vanishes when WG{t,x) is replaced by the exact solution of 
RG{uL,UR;to,XQ). Thus, the right hand side of (|2.8p represents the error due to the choice of 
approximate solution WG{t,x). 

Remark 2.2. 1. Condition (|2.7p ensures that the waves in Rc{uL,UR;to, xq) can not reach 
the lines = xo ± for t < to + s, so that the waves in the rectangle region D^j^^ ,^^-^ = 
[xq — r, Xo + r] x [tp, + s) do not interact with the waves outside Dfj-^^^^y 

2. In a different context, Liu [ID] derived earlier an estimate similar to (|2.7p . but for an 
approximation based on steady state solutions of the hyperbolic system and with initial data 
consisting of two steady state solutions of (2.1) (with / = f{u) and g = g{x,u)). 

3. Our formula (2.11) yields a possible generalization to the class of quasilinear systems 
(1.1) of the notion of (classical) Riemann solver introduced by Harten and Lax in [Hj . 

4. One can check similarly that Wg satisfies an entropy inequality associated with an entropy 
pair (when available). The error terms are completely similar to those found in (2.11). This 
will be used to show that the weak solution generated by the random choice method satisfies 
all the entropy inequalities. 

Proof. Without loss of generality, we can assume that {to,xo) = (0,0). Given a function 
9 with compact support in IR+ x K, we define m(i, x) := WgOiO + f{t, x, WG)dxd + g{t,x, VFg)^- 
From (|2.6p we have A(s, r; 9) = J^^ m{t, x)dxdt. Next, we decompose A(s, r; 9) as 



p 




(2.9) 



where 





and (if the i-wave, 1 < i < is a rarefaction wave) 




m(t, x) dxdt 
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We first compute in the region where classical Riemann solution Wc is a constant state. 
According to the form of Waitjx) in (|2.5p . it follows that 

WG{t,x)^u, + tq{0,0;u,) (2.10) 

for J e [a^ , CTj^i], i G {l,2,...,p— l}. By a simple calculation and the definition of q in ()1.4p . 
we have 

dtWa + d.J{t, X, Wg) - git, x, Wg) = g(0, 0; u,) - x, Wg) 

for j S {1,2, . . . ,p — 1}. By multiplication by the function 9 and then using integration by 
parts, we obtain 

A]{s,r;e) = / WGis,x)0{s,x)dx 



+ / (fit, <^i+it. WG{t, ai+,t)) ~ ar^^WGit, ar^itmt, ai+^t)dt 

(2.11) 

- / (/(t, a+i, WG{t, a+t)) - a+Wcit, a+t))e{t, a+t)dt 
Jo 

/ (q(0, 0; Wj) - q{t, x, Wg)) 0{t, x) dxdt. 

J aft 

By the property that q is Lipschitz continuous with respect to x and u on the compact set 
[0, s] X [— r, r] and the form of Wcit, x) in (|2.10p . the last term on the right hand side of (|2.1ip 
can be estimated by O(s^) ||^||co with the bound 0(1) depending on q. Therefore, equality 
(2.11) leads to 

A\{s,r-e) = / WG{s,x)e{s,x)dx 

+ ^ {fit. ^i+ii. WG{t, <7-+^t)) - a~^^WG{t, cjr^,t-))0(t, <7-^,t)dt ^2.12) 
{fit, a+t, Wait, a+t + 0)) - a+WGit, 4i)dt 



+ Oil)s^\\e\\co 

for i — 1,2, . . . ,p — 1. In the same fashion one can show that 

l-a-s M 

Alis,r;e) = / WGis,x)9is,x)dx ~ WGiO,x)0iO,x)dx 

J — T 'J — r 

+ ( ifit,a^t,WGit,a^t-))-a^WGit,a^t-))eit,a^t)dt 



(2.13) 



fit,-r,WGit,-r))eit,~r)dt 
Oil)s\s + r)\\e\\co+Oil)sr^\\e\\co, 



and 



Alis,r;d) = [ WGis,x)dis,x)dx ~ [ M^g(0, a;)6'(0, a;)da 

Jais Jo 

fit,r,WGit,r))0it,r)dt 



(2.14) 



if it, a+t, WGit, a+t+)) - a+WGit, cT+t+))e{t, a+t)dt 
+ Oil)s\s + r)\\9\\co + Oil)sr^\\e\\co. 
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Next, suppose that Wc{t,x) consists of an i-rarefaction wave in the region {{t,x)\j € 
[cTj^jCTj^]} for some i G 1, . . . ,p. It follows that Woitjx) in this region is of the form 

WGit,x) ^Wc{^) + tq{0,0;Wci^)) 

where Wc{j) is the i-rarefaction wave of the classical Riemann problem Rc{ul, ma; to^xo). By 

setting — J, Wcit.x) — WaitiC)^ the technique of change of variables {t,x) {t,S,), we 
obtain 

dtWc + dj{t, X, Wg) ~ git, x, Wg) 

= dtWG - ^d^WG + jdjit, t^, Wg) - g{t, t^, Wg) ^2.15) 
= \ {^{t, t^, Wg) - ^I) (/ + t|^(0, 0, Wc)) ■ ^ + g(0, 0,Wc) - git, t^, Wc) 

where I is the p x p identity matrix. Since Wc{£,) is a rarefaction wave for the system (|2.3p . 
this implies that 

Thus, by applying (|2.16p to (I2.15P we obtain 
dtWG + dj[t, X, Wg) - git, x, Wg) 

+ (^it,x,WG) ~ ^I)^{0,0,Wc) ■ ^ + q{0,0,Wc) - q{t,X,WG). 

Next, we multiply (j2.17p by 9{t, x) and integrate the equation over the region of z-rarefaction 
wave: t < s and j G [a~,a^]. Due to the Lipschitz continuity of ^ and the fact that 

f«' ^'^'-^ 1 remain bounded in [0, s] x [— r, r], the right hand side of p.l7p is bounded by 
0(l)s^(s + \ui — Ui-i\). Therefore, by p.l7p again, we deduce the estimate 

A2(s,r;6l)= / WGis, x)0{s, x)dx 



+ {fit, aft, WGit, a+t)) - a+WGit, att))0it, a+t)dt ^2.18) 

S 

if it, a-t. Wait, a-t)) - aiWcit, <y-t))eit, a-t)dt 



+ Oil)s\s+\uR-UL\)\\0\\co. 
Next, note that an j-shock wave satisfies the Rankine-Hugoniot condition 

/(0, 0, Ui) - (JiUi = /(O, 0, - aiUi^i, 
and this implies that the approximate solution WGit,x) satisfies 

[{fit, (T,t, Wait, (J^t+)) - a.Wcit, <y^t+))]e{t, a,t)dt 

- I [{fit, <y^-lt, WGit, a,_it-)) - a.^iWGit, a,^it-))]0it, a,_it)dt 
Jo 

= Oil)s^\un-UL\\\d\\co 



(2.19) 
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where the bound 0(1) depends on the Lipschitz constant of / and i°°-norm of q. Finally, by 
the estimates (EH), ([?J2)) - ([2J4)) and (|2lH)) - ((2l9)) . we obtain 

A(s,i;^) = ^Ai(s,t;0)+ ^ A2(s,t;0) 

/a~ s p— 1 r<7~j^iS 

WG{s,x)9{s,x)dx + ^ / WG{s,x)e{s,x)dx 
i=i •''^t^ 

+ ^ / WG{s,x)d{s,x)dx + / WG{s,x)d{s,x)dx 



VFG(0,a;)6'(0,a;)da; - / H/g(0, a;)0(O, a;)da; 



+ / fit,r,WG{t,r))eit,r)dt- f{t,-r,WG{t,^r))e{t,-r)dt 
Jo Jo 

+ 0{l){s^+r^){s + r + \uH-UL\) \\e\y, 
which leads to (|2.8p and completes the proof. □ 



3. Wave interaction estimates 

In this section we study the nonlinear interaction of waves issuing from two Riemann solu- 
tions and we derive estimates on the wave strengths. 

We emphasize that the generalized Riemann solution, nor the approximate solution Wcit, x) 
of the generalized Riemann problem RG{uL,UR',tQ,xo) is not self-similar. The solution does 
not consist of regions of constant value separated by straight lines. We thus should be careful 
in defining the wave strengths In fact, we still define here the wave strengths by using the 
underlying, classical Riemann solution Wc{t, x). We will see later that this strategy is accurate 
enough and that the discrepancy in total variation between WG{t, x) and Wc{t, x) on each time 
step is uniformly small (Cf. Section 4) when our Glimm scheme is applied to the problem (|l.ip . 
(|1.2p . The same observation applies to the potential of wave interaction to be introduced later. 

In the rest of the section, all waves are considered as waves from some classical Riemann 
problem unless specified otherwise. We say that an i-wave and a j-wave approach each other 
(or interact in the future) if either i > j, or else i — j and at least one of two waves is a 
shock wave. Suppose there are two solutions from different classical Riemann problems with 
strengths denoted by a = (a^, . . . , ap) and (3 — {f3i, . . . , f3p), then the wave interaction potential 
associated these two solutions is defined by 

^(a,/3) :=^|a,/3j|, (3.1) 

where the notation under the summation sign indicates an i-wave in one solution ap- 

proaching a j-wave in the other solution, and the summation is on all approaching waves; also 
ai or (3i is negative when i = j. In addition, given a {uL,Uf{;tQ, xq) € U x U x M+ x R, the 
wave strengths in Rg{ul, uh^Iq^xq) are denoted by e{uL,uii;to, xq). 
We first recall: 
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Lemma 3.1. (Glimm) 1) Given a (to, 2^0) in x M and ul, um, Ufj in U, we have 

\-1~{a + P)\=0{l)D{a,p) (3.2) 

where 

a ^ s{uL,UM;to,xo), P ^ e{uM,uii;to,XQ), j = e{uL,UR;tQ,xo). (3.3) 
2) Let vli vr he two constant states in lA, then 

D{j, S) = Dia, S) + DiP, d) + Oil)\d\D{a, l3), (3.4) 

and 

D{d, 7) = D{S, a) + D{6, p) + 0{l)\6\D{a, p) 

where a,P and 7 are given in (|3.3p . and S is given by 6 — e{vL,VR;to,xo). 

The following lemma describes the dependence of the wave strengths and potential D{-,-) 
with respect to their arguments. We introduce the following "local norm" of a given function 
(p{t,x,u) 

^tulTiv) = s^p[\V'(t,x,u)\ ;t e [ti,t2], xG[xi,X2], uGwj, (3.5) 
where the supremum is taken over any function u £ U and (t, x) G [ti, t2] x [xi, 2:2]. 

Lemma 3.2. 1) The wave strength e — (ei)i<i<p U U x M+ x ^ ^ W is a C'^ 
vector function of its arguments. Furthermore, for any {ul^ur), {u'j^,u'j^) in U x U and any 
(to, xq), (^0, x'q) in R+ X M, we have 

\a' - a| = 0(l)|a|(K - ul\ + \u'^, - ur\ + C^\t'o - to\ + C^lx'^ - xo\) 
+ Omu'j,-u'^)-iun-UL)\ 

where 

a ^ e{uL,UR;to,xo), a' ^ e{u'L,u'R;t'Q, x'q), (3.7) 
and the constants C'l and C2 are given by 

2) For given (ul,ur),{vl,vr),{u']^,u'j{), iv'^,v'j^) in U x U and {ti,xi),{t2,X2), {t'i,x[), 
(t'2, X2) in X R, we have 

D{a', P') = D{a,P) + 0{l)\a\\{v'R ~ v'^) - (vr - vl)\ 
+ Oil)\P\\iu'R~u'^)~{uR-UL)\ 

+ 0(l)|a||/3| (K -ul\ + \u'j^ - ur\ + K-vl\ + \v'r - vr.\) ^3 
+ 0{l)\a\\P\ {Cr\tL-tm\+C^'\xL~x^\} 

m=l,2 

+ Oil)\iu'j, ~ - (UR - ul)\ ■ Hv'r - v'l) - {VB. - vl)\ 

where 

a ^ e{uL,UR;ti,xi), P = e{vL,VR;t2, X2), {3 10) 

a' = e{u'L,u'j^;t[,x[), P' = e{v'^,v'^;t'^,x'2), 
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and the constants C™, C™ are defined by 



KZt-ij^), CT iV-->(^), m = 1,2. (3.11) 



Proof. The regularity of functions £i, i = 1, 2, . . . ,p, is a consequence of smoothness of the 
flux function / and the resuh of [15]. Moreover, the functions qiq^^ and q^q^^ are bounded if 



are bounded. 



dtdu ' dxdu 

To show p.6p . we note that ei{uL,uji]tQ^XQ) = when ur — u^. Then, by the regularity 
of £i, i — 1, 2, . . . ,p, we can express ei{uL, ur; to, xq), Si{u'i^, u'^r, t^, x'q) as 

/■I d^- 

ei{uL,uji;to,XQ) ^ / Tr-^(wL, (1 - t)wl + ru_R; to, a;o)dT • (ufl - u^), 
Jo (^""fl 

/■I 5^. 

e,«,ui^;to,a;o) = ^«,(1 -rX + ru^; ^q, a;^)dT • (u'^ - <). 
Applying the definition of {C* : j = 1, 2.} in (|3.8p and the norm in (|3.5p . we obtain 
ei(wL,Mfl;io'^o) £i{uL,UR;to,xo) 

(-— ^(w'i, (1 - t)u'^ + Tu'j^-tQ,Xo) - —-^{uL, (1 - t)ul + TUR;tQ,xo))dT ■ {ur ~ ul) 

OUr CIUr 

+ Q:^{'i^'LA'^-T)u'L+Tu'ji;t'o,Xo)dT ■ {(u'j^-u'l) ~ (ur-ul)) 

= 0(l){|u^ - ul| + \ur - ur\ + C°|<o -to\ + (7° l^o - a;o|}|u_R - ul\ 
+ Omu'j,~u'^)^{uR~UL)\, 

Therefore, by the observation of p.7p and the fact that 

Iwfl - ul| = 0(l)|£(uL,Ui?,;io, 2:0)1 = 0(l)|a|, 

we obtain p.6p . 

Next we derive p. 91) . By applying p.6p directly, we have 

a'^ = a, + 0(l)|a|{K - wl| + \u'r - ur\ + Cl\t[ - ii| + C^\x[ - 

+ O(l)|(li'«-0-(WK-"L)I, 

= /3j + 0(1)|/3|{|«1 ~vl\ + K- vr\ + - t2| + Ci\x'^ - X2\} 

+ Omv'^-v'^)-ivR-VL)\ 

for i, j = 1,2, . . . ,p where the constants {C™ : j,m — 1, 2.} are given in (|3.1ip and (|3.5p . We 
defined := {]< - ul| + lit^^ - + |i; - ii | + C]|a;i - xi]} and B {\v'l-vl\ + \v'r-vr\ + 
C^jtj — + C'||x2 — a^2|}- Then by multiplying two previous equations together and using the 
fact that A, B are of order 0(1) for {ul,ur), {u'i^,u'j^ GU xU, we obtain 

a'^P'j = a,/3, + 0(l)|a||/3|(A + B) + 0(l)|a||(z;^ - v'^) - (,;^ - vl)\ 
+ 0{l)\mn'R~u'^)~{uR-UL)\ 

+ 0{l)\{u'j^ - O - {uR - ul)\ ■ \{v'r - v'l) - {VR - vl)\, 

i,j = 1,2, ... ,p. Summing up previous equations for i, j — 1,2, ... ,p, we obtain (|3.9p . The 
proof is completed. □ 
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Using Lemmas 3.1 and 3.2, we obtain wave interaction estimates -which can be interpreted 
as a generalized version of ^lOj . 



Proposition 3.3. 1) Suppose that s, r are two positive numbers and (tQ,Xo) is in R_|_ x R. 
Also assume that ul, um, ur, ul + + ^lr are constant states in 14 and a, (3 and 7 

are the wave strengths of solutions of three classical Riemann problems Rc{ul, um', toj^o — r), 
Rc{uM,UB-,to,xo + r) and Rc{ul + fJ-L,UR + fj.B:,to + s,xa), i.e., 



(3.12) 



a = e{uL,UM;to,xo -r), /? = e{uM,UR;tn,xo + r), 

7 = S{UL + flL^UB. + fJ-Fi; to + S,Xo). 

Then we have 

|7| - \a\ + \P\+0{l)D{a,(3) 

+ 0(l)(|a| + \m\l^L\ + \i^r\ + C,s + C^r) (3.13) 
+ 0(l)|Mfl-/iL| 
where constants Ci and C'2 are defined by 

2) Let a, (3, 7 be the wave strengths as described in p.l2p . Also, for a given {vl, vu) inU x U 
and {ti, xi) in M_|_ x M, we define S = s{vl, vr] ti, xi). Then 



(3.15) 



(3.16) 



(3.17) 



-D(7, S) = D{a, S) + D{I3, 5) + 0{l)\5\D{a, [3) + 0{l)\5\\^lR - hl\ 
+ 0{l)\S\{\a\ + |/3|)(|ml| + \m\ + Cis + C^r), 

and 

D{5,-i) = D{5,a) + D{5,P) + 0{l)\5\D{a, (3) + 0{l)\6\W - 
+ 0(l)|<5|(|a| + \m\l^L\ + \iib\ + Cis + C^r) 

where constants Ci and C2 are given in (|3.14p . 

Proof. By the definition of 7 in (j3.12|) and Lemma 3.2 with u'j^ = ul + ij-l, 
u'p, — UR + /ifl,, t' = to + s, x'q = xq, wc obtain 

7 = e{uL,UR;to,XQ) + 0{l)\e{uL,uii]to,xo)\ (^\fiL \ + + Ci s^ 
+ 0{l)\nn,- fiL\ 

where constant Ci is given in (|3.14p . Similarly, by Lemma 3.2 we have 

£iuL,UM'ito,xo) = a + 0{l)C2\a\r, (3.18) 

eiuM,UR;to,xo) ^ P + 0{l)C2\f3\r. (3.19) 

On the other hand, Glimm's interaction estimates (|3.2p . (|3.3p lead to 

£{uL,UR;to,xo) ^ e{uL,UM;to,xo) + e{uM,UR;to,xo) 

+ 0{l)D{e{uL,UM;to,xo),e{uM,UR-to,xo)). 

Also, by (|3.9p - (|3.1ip with a' — e{uL,UM',to,xo) and 
/3' = £{um, ur] tQ,xo), we obtain 

D{£{uL,UM;to, xo), s{uM ,UR;tQ, Xq)) = D{a,l3) + 0{l)\a\\[3\C2r. (3.21) 
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Then, from ((3T7)l - ((33T1) it follows that 

|7| = \a\ + \p\+0{l)D{a,p)+0{l){\a\ + \p\ + \am)C2r 

+ 0{l)\e{uL,UR;to,xo)\{\^iL\ + + Cis) + 0{1)\^ir - , 
= \a\ + \p\+0{l)D{a,p) + 0{l){\a\ + \P\)C2r 

+ 0{l)\e{uL,UB:,to,xo)\{\^L\ + + Cis) + 0(l)|/ifl - 

Also, we see that estimates p.20p and (I3.2ip yield 

\£{uL,UR\t(i,xii)\ = \e{uL,UM]h,xo)\ + |£(um, Ui?,; to, 2:o)| 

+ 0{l)D{e{uL,UM\ to, xo),e{uM, ur; to, xo)) 

= \a\ + \P\ + 0{l)D{a, p) + 0(1) (|a| + Car 
+ 0(l)|a||/3|C2r 

= (|a| + + 0{l)C2r) + 0{l)D{a,(3) + 0(l)|a||/3|C2r, 

which in particular implies that 

\e{uL,UR-to,xo)\ - 0(1)(|«| + (3.23) 
Therefore, combining ([X^ with we obtain ((XT^ . 

Next we derive p.lSp . The proof of (|3.16p is similar, and is omitted. By the estimate (|3.4p 
we see that 



D{e{uL,UR:,to,xo),5) = D{e{uL,UM;to,Xo),S) + D{£{uM,UR;to,xo),5) 

+ 0{l)\5\D{e{uL,UM;to, xo),e{uM, ur; to, xo))- 

On the other hand, estimate (|3.9p yields 

D{l,5) = D{e{uL,UR;to,xo),S) + 0{1)\S\\^ir ~ 

+ 0{l)\e{uL,UR;to,xo)\\S\{\^j.L\ + \^J■R\ + Cis), 



(3.24) 



(3.25) 



and 



D{e{uL,UM; to, xo),S) = D{a, S) + 0{l)\a\\S\C2r, (3.26) 
D{eiuM,UR; to, xo),S) = D{(3, 5) + 0(1)|/3| \5\C2r, (3.27) 

D{e{uL,UM;to,xo),e{uM,UR;to,xo)) = D{a,f3) + 0(l)|a||/3|C2r. (3.28) 



Thus, by applying ([H?^ . (gSB-dSSHl) to ([gTM)) . we obtain the estimate ((XT5|) . The proof is 
completed. □ 

We just showed in Proposition 3.3 that Glimm's interaction estimates (Lemma 3.1) remain 
valid for the quasilinear hyperbolic system ()1.1|) up to certain error terms. The following 
immediate consequence of Proposition 3.3 will be the key to the forthcoming stability result. 

Corollary 3.4 Following the notations and assumptions in Proposition 3.3 and letting 
:= -sq{to + s, xo,ul), ^J-r := -sq{to + s, xo,ur) 
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in (3.14), we have 

|7| = \a\ + \p\+0{l)D{a,f3) 

+ 0{l){\a\ + |/3|)((Ci +C3 + C4) s + C2 r) , 

D((5,7) = D{S,a)+D{5,^)+Oil)\S\D{a,P) 

+ 0(l)|<5|(|a| + |/3|)((Ci + Cs + Ci) s + C2 , 

D{j,5) = D{a,S) + Dif3,S) + Oil)\S\D{a,(3) 

+ 0(l)|<5|(|a| + maC'i + C3 + Ci).s + C2r} 

where constants Ci, C2 are given in p.l4p and C3, C4 are given by 



(3.29) 



(3.30) 



(3.31) 



Cs N^:l\Mt, ^, N^ojo^^{^{t, x, u)). (3.32) 



Proof. By the observation of p.23p we obtain 

IfJ-R - fJ-L\ = s\q{ta + s,XQ,uii) - q[tn + s,xq,ul)\ 

dq , _, . . 

= s—(to + s,xo,u) ■ \uii- ul\ (3 33) 

= 0{l)C4,s\e{uL,UR;to,xo)\ 
= 0{l){\a\ + \(3\)C4S 

where u ^ U and C4 is given in (I3.32p . Therefore, by combining (|3.33p with the resuh of 
Proposition 3.3, we obtain (|3.29p - (|3.3ip . The proof is completed. □ 



4. Stability of the generalized Glimm method 

We are in position to introduce our version of Ghmm scheme for the approximation of the 
quasihnear system p.ip . Then we rely on the wave interaction estimates in Section 3 and prove 
a stability result. 

The approximate solution to the Cauchy problem (|l.ip . (|1.2p is defined as follows. Given 
two positive constants s and r satisfying the C-F-L condition (|2.7p . we introduce the constant 

A, := -. (4.1) 

s 

Let also a — {a^. : oj. G (—1, 1),A; G N} be an equidistributed sequence. We divide the {t,x) 
plane into 

tk^ks,Xh^hr, fc = 0,1,2,..., /igZ. (4.2) 

Next, we construct an approximate solution Ur{t, x) of the problem (|l.ip . (|1.2p in the following 
way. First, the initial data uo{x) is approximated by a piecewise constant function 

Ur{0, x) = uo{hr), X e[{h- l)r, (h + l)r), h is odd. (4.3) 
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Then, within domain < t < s, we construct an approximate solution Wcit^x) for each 
generalized Riemann problem with initial data Ur{0, x) to obtain Ur{t, x) in region x); < 
t < s}. If Ur{t, x) has been constructed for t < ks, A; G N, we set 

Ur{ks,x) :— Ur{ks—, {h + ak)r) (4-4) 

for a; G [{h—l)r,{h + l)r), k + hisodd. Again, we solve the generalized Riemann problems with 
initial data Ur{ks, x) given in (14. 4p to construct Ur{t, x) within region {{t,x);ks < t < {k + l)s}. 
Following the process (14. 3p . (j4.4p consecutively, we then construct our approximate solution 
Ur{t,x) of (|l.ip . (|1.2|) . In other words, the approximate solution to the problem (|l.ip . (|1.2p 
generated by the generalized Glimm scheme is given by 

Ur{t, x) — Wait, x; Ur{ks, [h — l)r), Ur{ks, {h + l)r); ks, hr) (4-5) 

for {t, x) G [ks, [k + l),s) x [{h — l)r, [h + l)r), fc + is even. 

Next we study the stability of Ur(t,x) in L°° and BV^ norms. This requires the description 
of mesh points, mesh curves and immediate successors beforehand. Recall that the values of 
Ur{t,x) on t = ks are determined by the values of Uj.{t,x) at points {{ks—, {h + afc)r); h G 
Z, A: + ft, is odd}, we call these points {{ks, {h + ak)r) : k — 0,1,2, - ■ ■ , h e Z, k + h is odd} the 
mesh points of approximate solution Ur{t, x). We obtain a set of diamond regions by connecting 
all mesh points with segments. An unbounded piecewise linear curve / is called a mesh curve 
if / lies on the boundaries of those diamond regions. Suppose / is a mesh curve, then / divides 
the {t, x) plane into /"'" and /~ regions, such that /~ contains i = 0. We say two mesh curves 
Ii > I2 { Ii is a successor of I2) if every point of Ii is either on I2 or contained in I2 ■ And, Ii 
is an immediate successor of I2 if Ii > I2 and every mesh point of Ii except one is on I2 ■ Note 
that the difference between Ii and I2 is determined by a diamond region if one is an immediate 
successor of the other. 

Next, to simplify the notations, we set U]~^h '■= Ur{ks,hr) when k + h is odd. By the 
observation of (|2.4p and (|4.5p , we have 

UkM = Uk,h + s q{{k - l)s, hr, iik,h), k + h is odd, 

where Uk,h is the value of Rc{uk-i,h-i,Uk-i,h+i; {k - l)s, hr) at {ks-, {h + ak)r), i.e., 

r 

Uk.h = Wc{ak-\Uk-i.h-iiUk-i.h+ii [k - \)s, hr) 
s 

with the function Wc given in Section 2. Next, given a pair (feg, h^), kg + ho is even, we note 
that the (i, a;)— plan consists of the diamond regions Tkg,ho with center {kos,hor) and vertices 
(mesh points) 

S := ((fco - l)s, {ho + ak„-i)r), W :== {kos, {ho - 1 + afejr), 
E := {kos, {ho + l + ak„)r), N := {{ko + l)s, {ho + ak„+i)r) 

(see Figure 4.1). We set 

Us ■— Uko-lMm Uw '■= Uk„,ho~l: Ue ■= Uk„,hg + 1, "TV := Ukg + lMo, (4-7) 

and 

ILS '■— Ukg-lMo, Uw '■= Uka,ho-l, UE '■= Ukg.ho + l, UN ■= Uko + lMo- (4-8) 

Note that uw and ue are the states in Rc{{ko — l)s, {ho — l)r) and i?G((fco — l)s, (^o + 1)^) 
respectively, i.e., 

Uw = uw + s q{{ko - l)s, {ho - l)r, Uw), 
ue = UE + s q{{ko - l)s, {ho + l)r, He)- 



16 



t = (irg + l)s ■ 



N 



t = k^s 











1 ^ X 













Figure 4.1 : Diamond region Tko.ho 

Now we define the strengths of waves in Ur{t, x). However, the set up for the waves strengths 
of Urit^x) becomes crucial due to the lack of self-similarity of approximate solution Wcit^x), 
the strengths of waves in Wcitjx) can not be defined in the traditional way as described in 
[15j . To overcome the difficulty, we first solve the associated classical Riemann problems with 
the initial data {ur{ks—, [h + ak)r); x G [{h — l)r, {h + l)r). A: + ft, is odd} (see (|4.4p ) within 
each time step. So we construct a new function Ur{t, x) defined on R+ x M. Then we define the 
strengths of approximate waves in Ur(t, x) based on classical waves in Uj.{t, x). More precisely, 
given a wave in Ur{t,x), there exist two corresponding constant states Ui 

and a classical Riemann wave {ui-i,Ui) with strength e{ui-i,Ui) in Ur{t,x), then the strength 
of {ui-i{t),Ui(t)) is defined as 



Next, we show that, under the condition that the L^(I 



-norms of q and ^ are small 



the sum of strengths for waves in Ur{t, x) crossing mesh curve J can be regarded as an equivalent 
norm for the total variation of a;) on J. By the fact that the term \e(ui-i,Ui)\ is equivalent 
to the total variation of {ui-i,Ui) for any classical Riemann wave it is equivalent to 

show that the total variation of Ur{t, x) on J is equivalent to the total variation of Mr-(i, x) on 
J. To show this, let Jk be a mesh curve lying within fc-th time level {{t, x); ks < t < {k + l)s}, 
and let TV{ur{t, x), Jk), TV{ur{t, x), Jk) denote the total variations of Ur{t, x), Urit, x) on Jk 
respectively. Suppose there is a wave {ui-i{t), Ui{t)) in Ur{t, x), issued from {ks, ir) and crosses 
Jk, also is the corresponding classical Riemann wave of (ui-i(t),Ui{t)) (so (ui-i,Ui) 

is also issued from {ks,ir) and crosses Jk)- If (ui-i,Ui) is a shock wave, then by (|2.4p we can 
easily obtain that 

\TV{{u,^i{t),u^{t)); Jk) - TV{{u,^i,u,); Jk)\ 

< s\ — {ks,ir, Ui)\ TV{{ui_i,Ui); Jk) + s (|g(to, xq, + \q{to,xo,Ui)\) , 

where Ui € U and TV((ui^i{t),Ui{t)); Jk), TV{{ui^i,Ui); Jk) denote the total variations of 
(ui^i{t),Ui(t)), (ui-i,Ui) crossing Jk- Similarly, if {ui^i,Ui) = Ui{C) is a rarefaction wave with 
C G [^iiCa], then we obtain 



\TV{iu^-lit),U^it))■, Jk) - rF((w,_i, U,); Jfe)l 

dq 

< s\ — {ks, ir, Ui{£,)\TV{{ui-i,Ui); Jk) + s{\q{to, Xq, Ui-i)\ 
ou 



\q{tQ,xo,Ui)\) 



for some ^ £ [^1,^2] and Ui{^) G U . Summing up the previous inequalities with respect to the 
waves crossing Jk we obtain 



\TViuriJk))-TViMJk))\ 
dq II 



<Ois) 



TV{Ur{Jk)) + 0{s) ||g||Li(R^xR) 
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for any mesh curve Jfe, and this is enough to imply that the total variations of Ur{t,x) and 
Ur{t,x) on any mesh curve Jk are equivalent when ||(?||li(r^xR) and || |^ ||li(r+xR) are small, 
we then show the statement. 

We note that the waves entering each diamond region may come from two generalized 
Riemann solutions, we certainly need to know the constant states of corresponding classical 
Riemann solutions at the left and right vertices of diamond region to calculate those wave 
strengths separately. We proceed as follows. 

First, using the notations in (|4.7p . (|4.8|) . we define the strength of the waves entering the 
diamond region Vkg^ho, + ho is even, by 

£*{Tko,ho) ■= \£{uw,us; {ko - l)s, {ho - l)r)| 
+ \s{us,ue; [ko - l)s, [ho + l)r)\ 

and the strength of the waves leaving T^g^ho by 

£*(rko,ho) \s{uw,UN;kos,hor)\ + \e{uN,UE;kQS,hor)\. (4.9) 

Since uat is a constant state in Wc{uw, UE',kQS, hor), we can write 

£*{^koJio) \£[uw ,UE;kos,hor)\. (4-10) 

Next, for ko + ho is even, we let Q{Tko.ho) denote the potential of waves interaction in the 
diamond Tko.ho^ i-S-, 

Qi^koMo) D{e{uw,us, {ko - l)s; {ho - l)r),e{us,UE, (ko - l)s, (ho + l)r)) 

where D(-, ■) is defined in (|3.ip . Given a mesh curve J, we note that there are two types of waves 
crossing J. The first kind of waves are (uk,h-i,Uk-i,h), k + h — even (waves of type I), the 
second type of waves are {uk-ijn Uk,h+i), k + h = even (waves of type II). More precisely, waves 
of type I are either of the form {uk,h-i,'Uk-i,h) entering Tk,h (left in-coming waves ofTk^h), or 
{uk+i,h, Uk,h+i) leaving Tk^h (right out-going waves of Tk.h)- Waves of type II are either of the 
form {uk-i,h,Uk,h+i) entering TkM (right in-coming waves ofVk^h), or (uk,h-i,Uk+iM) leaving 
Tfe.ft (left out-going waves of Tk.h)i see Figures 4.2 (a), (b). Next we define the linear functional 
L{J) for the waves in Ur{t, x) crossing mesh curve J by 

L(J) X! \£('>^k,h~i,Uk-i.h;(k - l)s,(h - l)r)\ 

type I , , 

^ (4.11) 

-I- 2^ \e(uk-i,h,Uk,h+i;{k-l)s,(h+l)r)\. 

type II 

From previous analysis, we see that functional L{J) is equivalent to the total variation of u,.(i, x) 
crossing mesh curve J. Next we define the quadratic functional Q(J) of Ur{t, x) by 

Q(J) ^ i?(a,/3) (4.12) 

(",/3) 

where the notation {a, [3) under summation sign denotes a pair of waves a, (3 crossing J and 
approach, and D{a,l3) is given in (|3.ip . Furthermore, we define the Glimm functional F(J) 
of Ur(t,x) for mesh curve J by 

F{J) := L(J) + KQ(J). (4.13) 

Our goal is to show that functional F remains uniformly bounded on all mesh curves provided 
that constant K in (|4.13p is sufficiently large, and this leads to the result that functional L 
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can be bounded by a constant times the total variation of initial data uo{x). To show this, the 
first step is to estimate the possible changing amount of L and Q when waves pass through one 
mesh curve and into an immediate successor. The estimates of changing amounts of L and Q 
are stated as follows. 



Figure 4.2(a) : Waves of type I crossing mesh curve J 



Figure 4.2(b) : Waves of type II crossing mesh curve J 



Proposition 4.1. Given two mesh curves Ji and J2 such that J2 is an immediate successor 
of Ji, let ^ko,ho denote the diamond region bounded by Ji and J2. Then Junctionals L and Q 
satisfy 

L{.h) - L{.h) = 0{\){Q{Tk,M,) + e*(rfeo,/.o)(C? + + Cl + C°)s}, (4.14) 



Q{J2) - Q[Ji) = - Q{TkoM) + 0{l)L{J,)Q{Tk„M) 

+ 0{\)L{.h)e,{Vk,M){C'i + KCl + Cl + Cl)s 



where constants A, is defined in (|4.ip and ,1 < j < A, are given by 



C! N^^;XMMt^-^-))^ Cl := N^^orMr^^^j2{t,x,u)) (4.17) 
where N is defined in (|3.5p . Note that {Cj*; 1 < J < 4} depend on ho,kQ. 

Proof. Let us,uw,ue,un be the constant states described in (|4.6p - (|4.7p . we first derive 
(|4T4|) . By the definitions of e, and e* in gj]), (|4T0l) and L in ((4TT|) . we find 

L(J2) - L{Ji) = \e{uw,UN;kos,hQr)\ + \e{uN,UE;kos,hor)\ 

- \e{uw, us; (fco - l)s, {ho - l)r)| 

- \s{us, ue; (fco - l)s, {ho + l)r)\ 
= £*(Afco,/iJ - e^{AkoMo)- 

Next, by applying the definition of A, in (|4Tl) and the estimates (|3?29l) . (|3?32| to (|4T8ll with 
the choice of ul = uw, um — us, ur — He, /^l — uw — uw, fJ^R. — ue — ue, to = (fco — l)s 
and xo — hor, we obtain 

e*{rk„,ho) = e.{TkoM) + 0{l)Q{rk„M) + 0(l)£*(rfe„,,j((C? + Cl + Cl)s + C°r) , 
and this gives (|4.14p . 
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To prove (|4.15p . we define several notations for the rest of the section. First, given {k, h), k + 
h= even, we let vector e^.i /j_i/2 denote the strength of waves issued from ((fc — l)s, {h — l)r) 
entering Tk.hy and let vector efe-i,/i+i/2 denote the strength of waves issued from ((/c — l)s, (h + 
l)r) entering Tk^h- More precisely, the vector £fe_i.;t_i/2 measures the strength of waves of 
type I entering Tk^h and ek-i^h+1/2 measures the strength of waves of type II entering Tk^h- 
Next, given a mesh curve J, let J[h-i.h] {Jlh.h+i] respectively) denote the segment of J in 
M+ X [{h — l)r,hr] (R+ x [hr,{h + l)r]). Then we define vectors ej^h-1/2, £j,/i+i/2 as the 
strengths of waves crossing J[{h-i).h]j Jlh,h+i] respectively. We will drop the sign J in £j^h-i/2 
and £j^h+i/2 when J is specified. We also set 

£W,S ■— £ho-l/2, £S,E ■— £/io + l/2i 

ew,N := £iuw,UN; kos, hor), sn.e := £{un, ue; kos, hor). 
Since J2 is an immediate successor of Ji , the diamond region bounded by Ji , J2 can be specified 
as r(fco)^o) with center (/cps, hor), and Ji, J2 coincide outside r(A:o,^o)- We will also drop 
the signs Ji, J2 without confusion. From the definition of Q in (|4.12p . we have 

Q(J2) - g( Ji) 

= ^ (^D{eh-i/2^ £w,n) + D{eh_i/2, £n,e) - £'(£/i_i/2, £w,s) ~ D{efi_if2,£s,E)^ 

h<ha 

+ X! {p{ew,N, £h-i/2) + D{eN,E, e/i-1/2) - D{£w,s, £11-1/2) - D{es,E, £^-1/2)) 

h>ho + l 

+ D{ew,N, £n,e) — D{ew,s,£s,E)- 
From (j3.1|) we see that 

D{ew,N, £n,e) = 0. 
Also, for any ft, e Z we observe that 

D{£h-i/2, £w,n) + D{£h-i/2, £n,e) = 

for h < ho, and 

D{£w,N, £11-1/2) + D{eN.Ei £h-l/2) = 
for h> ho + 1. Thus, by (|lTl])-(|12Il) we obtain 
Q{J2)-Q{Ji) 

= ^ (^D{eh-i/2,£w,E) - D{eh-i/2,£w,s) 

h<ha 

+ (D{ew.E ■,£h-i/2) - D{ew,s,£h- 

h>ho + l 

- D{ew,s,£s,E)- 

Finally, applying p.30[) and p.3ip to (|4.22p and using the fact that D{ew,s,£s,E) = (3(rfco,/io)' 
we obtain 

Q{J2) - Q{Ji) 

= -D{ew,s,£s,E) + Y (o{l)\£h-i/2\D{ew,s,£s,E) 

h^ho,ho + l 

+ 0(l)|£/.-i/2|(|ew.s| + \£sM)iiC°i + C° + C"^)s + C°r)) 
= -Q(rfc„,,.J +0(l)i(Ji)Q(rfc„,,J + 0(l)L(Ji)e,(rfc„,,.J((C? + C3VQ")s + C°r), 
which leads to (|4.15p . This completes the proof. □ 



D{£h-1/2,£W,E) 

D{ew,E, £1-1-1/2) 



(4.19) 
(4.20) 
(4.21) 



- D{£h-1/2,£S,E) 

-1/2) - D{es,E, £h-i/2, 



(4.22) 
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Before stating a crucial technical lemma, let us introduce a notation about mesh curves. 
We say that a mesh curve J is of the type {ko, ko + 1) if all the mesh points on J have the form 
of {{ks, {h + ak)r) : k — ko, ko + 1}. 



Lemma 4.2. Given a positive integer ko, let Ji and J2 he two mesh curves of type (feg — 1, ko) 
and {ko,ko + 1) respectively. We assume that there exists a positive constant such that 

L{Ji)<M,. (4.23) 

// is sufficiently small and the constant K in (|4.13p is sufficiently large, then the functional 
F satisfies the following inequality 

F{J2) < F{J,) + 0{l)s ^ e,{TkoM)iCi + KC^ + + C°) (4.24) 

hoez 

where the bound 0{1) depends on M^, and K, and the constants Cj := Cj{ho,ko), 1 < J < 4, 
in ([4T61) . (|4T7l) depend on ho £ Z. 



Proof Given ho 6 Z, we multiply (|4.15p by constant K in (|4.13p and add it to (|4.14p . Then 
by the assumption that Ji and J2 are two mesh curves of type {ko — 1, fco) and {ko, ko + 1), we 
obtain 

F(J2)-F(Ji)= ~Kj2Qi^ko,ko) + 0{l)[l + KL{J,)]i 5]Q(r,„,J 
+ J2 ^*i^'^oM)iC'l + KC^2 + + Cl)s\. 

Next, by the observation that X^/mez Q(^*:o.'io) = Q{Ji)i the equation above implies that 
F{J2) - F{Ji) = -KQ{Ji) + 0(1)(1 + KL{.h))Q{.h) 

+ 0{l)s e*{^koMo)iC°i + KC^2 + C3" + Cl) 

=.Q(Ji){i^[0(l)i(Ji)-l] + 0(l)} 

+ 0{\)s J2 £*{rkoM){C'^ + A.CO + C3" + Ci) 

hoez 

<Q{Ji){K[0{l)AU-l] + 0{l)} 
+ 0{l)s J2 e.{Tko,ho){C'l + KC^2 + C3" + Cl). 

haeZ 

The last inequality is an application of (|4.23p . We see that the term K[0{\)Mf. ~ 1] + 0(1) is 
negative, if is sufficiently small and K is sufficiently large. Thus, (|4.24p holds for such 
and K. This completes the proof. □ 

We now establish the stability of generalized Glimm method, which is the main result of 
this section. We denote by TV{-) the total variation of a function. 



Theorem 4.3 Fix a constant state and assume that the initial data uq — uo{x) is a 
function of bounded variation such that 

\\uo — u^,\\l°^ and TV{uo) are sufRciently small. (4-25) 
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Assume also that the mappings A{t,x,u) :~ ^^{t,x,u) and q{t,x,u) in (|1.4p are smooth and 
such that 

d^A d^A dq 
the L^(M+ X M) norm of „ „ , A*———, q, — are sufRciently small. (4.26) 

otau axon an 

Then, the approximate solutions Ur{t,x) are bounded uniformly in the L°° and BV norms: 

\\ur - '«*||l=°(r+xR) < 0(1) {\\uq -u*||loo(r) + TV{uq) + , (4.27) 



TV{ur{t, •)) < 0(1) {tV{uo) + C) , 



where 



C := 



d'^A 



dtdu 



LifR+xI 



+ A* 



d'^A 



dxdu 



+ lkllLi(R+xR) + 



dq 



du 



LH 



(4.28) 



(4.29) 



Furthermore, the function Ur{t,x) is Lipschitz continuous in time, i.e., for ti,t2 > 0, 

\ur{ti,x) - Ur{t2,x)\dx < 0(l)(|t2 ~ h\ + s) {TV{uo) + C) . (4.30) 



Proof. We apply an induction argument based on Lemma 4.2 to show that the approximate 
solution Ur{t,x) is uniformly bounded in L°° and total variation. First, we show that the 
condition (|4.23p in Lemma 4.2 holds under the assumptions (|4.25p . (|4.26p . By induction, given 
ko G N, we let Jko-1/2 denote the mesh curve of type (fco — 1, kg). For fco = I7 we see that 



F{Ji/2) < 0(1) {TV{uo)+K[TViuo)]' 



(4.31) 



This means that there exists a positive constant M,, as described in (|4.23p . such that F(Ji/2) < 
M,, and in particular, L(Ji/2) < if TV{uo) is sufficiently small. Next, suppose that 



L{Jk+i/2) < M* for fc - 0, 1, . . . , fco - 1. 



(4.32) 



We intend to show that (|4.32p still holds for k — fco- Since Jko-1/2 is a mesh curve of type 
(fco — 1, fco), this implies that Jka+1/2 is a mesh curve of type (fco, fco + 1) so that Lemma 4.2 
can be applied. Therefore we obtain 

F{Jko+i/2) < F{Jk„-i/2) + 0{l)s J2 e.iTk„M)iC°i + A*C° + C° + C°) 



Then, by 



this leads to 



< F{J,/2) + 0(l)s^ e*(rfc./.o)(Oi" + A,C° + 0° + C°). 
fc=i hoez 



e,iTk,ho) = L{.Jk-i/2), keN, 



kg 

FiJko+1/2) < F{J,/2) + Oil) J2 sL{Jk-i/2) sup (C° + X.C"2 + C3" + Q"). 



(4.33) 
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Next, by Km - Km we find 



F{Jko+i/2) < + KTV{uo)) TV{uo) 

+ 0(1)M, sup (C? + A*C2° + + C2)s. 

From the definitions of and the assumption that the constant C in (|4.29p is finite, we see 
that 

oo 

lim y sup (Ci° + KC^ + C3" + C2)s = C. (4.35) 

Therefore, from (|4.34|) and (|4.35p we obtain the inequality 

^^(^■0+1/2) < 0{1){{1 + KTV{uo))TV{uo) + M.C}, (4.36) 

and in particular, 

L{Jk,+i/2)<0{l){{l + KTV{uo))TV{uo) + M,C}. (4.37) 

We note that the functional L in (|4.37p only depends on the constants M,, C and the total 
variation of uq, thus it enables us to choose TV{uo) and C sufficiently small such that 0(1)(1 + 
KTV{uo))TV{uo) < ^ and 0(1) CM., < ^ and this implies that 

LiJko+1/2) < M,. 

Therefore (I4.32p holds for k ~ ko, we just showed that L{ko + 1/2) has uniform bound for all 
fco G N, which implies that functional L of Ur(t,x) has global bound. Since L is a functional 
equivalent to the total variation of Ur{t.,x), we prove that the total variation of Ur{t,x) has 
an uniform bound for all i > and all finite r > 0, so as well the L°° norm of Ur{t,x). 
To prove (|4.28p . we apply (j^i]) . (|4?5|) to Ur{t,x) and we use the fact that TV{urikos, ■)) = 
0(1) i^( Jfco+1/2) to ([Qe]) . then (jiTM]) is established. For the proof of (14271) and ^M), we 
follow the lines of proof in [10]. The proof is completed. □ 

We note that if A,-^^, q and |f in (|4?26l) belong to L°°, then inequalities (f427|) . 

(|4.28p and ()4.30p remain valid in a finite interval [0, T] with T sufficiently small. 



5. Convergence of the generalized Glimm method 

In Section 4 we established the BV stability of the scheme together with a time continuity 
property. By Helly's theorem , there exists a subsequence of approximate solutions, still denoted 
by {ur{t,x)} and converging strongly in Lj^^ to a limit function u = u{t,x). Moreover, by the 
estimates ()4.23p . (j4.24p . the function u is uniformly bounded and is of bounded variation in x. 
We now prove that the limit u is indeed an entropy solution of the Cauchy problem. The proof 
relies on the error estimate derived in Section 2. 

Theorem 5.1 Suppose that the initial data uo{x) is sufficiently close to a constant state in L°° 
and BV , and that the norms of g,j.^^ , , q, and ^ are sufficiently small in x M. Let 
{itr(i, x) : r > 0} 6e the sequence of approximate solutions constructed by the generalized Glimm 
scheme (4-3)-(4.5). Then, for any equidistributed sequence {ak}keKj there exists a subsequence 
of {ur{t,x)} converging in L\^^ to a function u = u{t,x) which is an entropy solution of the 
Cauchy problem (1-1), (1-2). 
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Remark 5.2. Assume that U is a convex subset of R^, we say that {U, F), U : U ^ MP 
and F : M.^ x M x W ^ R, is an entropy pair of the system (|l.ip if [/ is a convex function on U 
and 

dF DU df ^ ^ 

— = on ]R+ X M X W. 

ou Du du 

Furthermore, a function u : x M — > is called an entropy solution of p.ip if u = u{t, x) is 
a weak solution of (|l.ip satisfying 

dtU{u) + c),(f^(i, X, u)) < ^{u){g{t, X, u) ~ {d^m, X, u)} + x, u) (5.1) 

in the sense of distributions, for every entropy pair {U, F). 



Proof. The proof is based on the result of Proposition 2.1. Let {ur{t,x)} denote a sequence 
of approximate solutions constructed by generalized Glimm scheme (I4.3l) - (|4.5p . Then, by the 
stability result and Helly's theorem, there exists a subsequence of {ur{t, x)} converging almost 
everywhere to a function u G Lj^^ with bounded total variation. Given any test-function 
9 : M+ x K — > R with we define the residual of 

TZ{ur,9) :— i i {urdtO + f{t,x,Ur)dxO + g{t,x,Ur)9}dx dt + / uo{x)9{0, x)dx 

JR+ JtS. J-oo 

Note that u is a weak solution to (|l.ip . (|1.2p if and only if TZ{u, 6) — for any test-function 9. 
By Lebesgue's theorem, we see that 

\n{ur,9)-n{u,9)\-^0 asr^O. 

Thus, to show that u is a weak solution of (jl.ip . (|1.2p . it is equivalent to show that TZ{ur{t, x), 9) 
tends to zero as r vanishes. To show this, we first let y'^"'''" denote a characteristic function 
having the same support as the test-function 9. Then, by construction of Ur{t, x) and by ([2 
we can write 

TZ{Ur,9) 

oo Ako + l)s Aho + l)r 

= / / (urdt9 + f{t,x,Ur)dx9 + g{t,x,Ur)9jdxdt 

ko=0 i^o+i^o J{ho~l)r 



= E E 0{l){s' + r^){s + r+\uk,.M+i-UkoM-i\)x'::pp\o) 

ko=0 ''o + '^o 

°° / i-(ho + l)r 

V V / Ur{{ko + l)s-,x)9{{ko + l)s,x)dx 

ko=0^o+^o V^(/.o-l)r (5.2) 
(ho + l)r \ poo 

Ur{kos+,x)9{kos,x)dx \ + / uo{x)9{0,x)dx 



{ho-l)r 

{ko + l)s 



EE/ f{t,{ho + l)r,Ur{t,{ho + l)r-))9{t,{ho + l)r)dt 

{k„ + l)s \ 

fit, {ho - l)r, ur{t, {ho - l)r+))9{t, {ho ~ l)r)dt . 

kos I 
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Let ^i{r), fl2{r) and flsir) denote the terms on the right hand side of (|5.2[) . respectively. We 
first estimate ^i{r). By a direct calculation and (|4.2[) . (|4.28p . we obtain 



n^{r) =0(l)r + Oa){s^ + r'){\ukMi - ^koM-iDxtpp" (e) 



kn=0 '»0+'=0 



<0(f)r+ 0{l){s'+r^){T.V.{uo{x)} + C)x'::,^; 
<0(l)r. 



(5.3) 



supp (O) 

kaEN 



Next we calculate n^{r). By the property of the Lipschitz continuity of /, g and 
we obtain 

= OWE E / K(i,(/io + l)r+)-zi.(i,(/io-l)r-)|-(x^:;f;(,))dt 

= 0(1) E E ( / ^l9(^os, (/lo + 2)r,fife„,,,„+i) - g(fcos,/ior,Ufe„ji„+i)| • {xTupp ^g))dt 



= o(i)Ee/ 



(fco + l)s 



It follows that 

ns{r)=Oil)r. (5.4) 
It remains to estimate (?') • It is a standard matter to check that 

p(/io + l)r 



f^2(r-) E / K](fcos,x)0(fc, 



'os, x)(ix 



+ 00 

{ur{0, x) ~ uo{x))6{0, x)dx 



where [ur\{kQS,x) := Ur{kos+,x) — Ur{kQS—, x). We let J{{ak},r,9) denote the term 



°° Aho + l)r 

/ [u,.](fcos, a;)6'(fcos, a;)(ia;. 

7. _i J(hn — l)r 



By the construction of Ur{t,x) in (|4.3p . we see that the term J^^{ur{0,x) — uo{x))9{0, x)dx 
on the right hand side of (|5.5p vanishes as r tends to zero. In addition, by a result of Liu 
[TO] we obtain that, for any equidistributed sequence {a^jfegN, J({afe},r, 0) tends to zero as r 
approaches to zero. This implies that 

r22(r)^0 asr^O (5.5) 

for every equidistributed sequence {0^}^^^. We refer the reader to [TO] for the details of the 
estimate of 5l2(»')- Finally, by (|5.3p . (|5.4p and (|5.5p . we obtain 

7^(wr,6') ^ inL^ as r ^ 0, 
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which means that the hmit function u satisfies TZ{u, 9) = 0. Therefore, u is a weak solution of 
the Cauchy problem (fLTj) . 

To prove that u is an entropy solution satisfying the entropy inequality (15. ip , it is equivalent 
to show that, for any entropy pair (U, F) and test-function 9 > 0, the function u satisfies 

U{u)9t + F{t,x,u)e^ + P{t,x,u)9dxdt + / U{uoix)) 9{x,0) dx > 0, (5.6) 



with 

P(t, X, u) := ^ • {g - |^)(^, X, u) + {dxF){t, x, u). 

We note that the result of Proposition 2.1 can be applied to show that u{t,x) satisfies (|5.6|) 
for any entropy pair {U,F). In turn, this implies that u is an entropy solution of the Cauchy 
problem (|l.ip . (|1.2p . and the proof of Theorem 5.1 is completed. 

□ 
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